Front explosion in a periodically forced surface reaction 
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Resonantly-forced oscillatory reaction-diffusion systems can exhibit fronts with complicated in- 
terfacial structure separating phase-locked homogeneous states. For values of the forcing amplitude 
below a critical value the front "explodes" and the width of the interfacial zone grows without bound. 
Such front explosion phenomena are investigated for a realistic model of catalytic CO oxidation on 
a Pt(110) surface in the 2:1 and 3:1 resonantly- forced regimes. In the 2:1 regime, the fronts are 
stationary and the front explosion leads to a defect-mediated turbulent state. In the 3:1 resonantly- 
forced system, the fronts propagate. The front velocity tends to zero as the front explosion point 
is reached and the final asymptotic state is a 2:1 resonantly- locked labyrinthine pattern. The front 
dynamics described here should be observable in experiment since the model has been shown to 
capture essential features of the CO oxidation reaction. 

PACS numbers: 82.40.Bj,82.40.Np,05.45. Jn 
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INTRODUCTION 



Catalytic surface reactions are a large class of chemi- 
cal reactions with important applications. They are of- 
ten accompanied by oscillations, and a rich variety of 
wave patterns, as well as chemical turbulence, can be 
observed in such systems 0, Q- An extensively stud- 
ied surface reaction is CO oxidation on single crystals of 
platinum. When the reaction takes place on the crys- 
tallographic plane Pt(110), a structural phase transition 
in the top layer of the substrate is coupled to the reac- 
tion, making oscillations and excitability possible. These 
kinetic regimes are well reproduced by the theoretical 
model of Krischer, Eiswirth and Ertl |3( where reaction- 
induced surface reconstruction is taken into account. The 
bifurcation analysis of this model shows that it has a su- 
percritical Hopf bifurcation |3|. Moreover, in a region 
of parameter space the uniform oscillations in this sys- 
tem may be unstable with respect to the modulational 
Benjamin- Feir instability, leading to the spontaneous de- 
velopment of turbulence |J. Therefore, in the vicinity 
of the Hopf bifurcation, the CO oxidation reaction ap- 
proximately obeys the complex Ginzburg-Landau equa- 
tion (CGLE) 5]. Experimentally, it may be difficult to 
maintain the reaction so close to the bifurcation line that 
the oscillations have a harmonic character and their am- 
plitude is small. Nonetheless, many qualitative features 
of the wave patterns typically found in the CGLE are 
also observed in the CO oxidation reaction on Pt(110) 
under oscillatory conditions. In particular, the sponta- 
neous development of amplitude turbulence through the 
nucleation of spiral wave pairs has been found in this 
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system 0, S 0- Thus, the CO oxidation reaction on 
Pt(110) provides an opportunity to experimentally inves- 
tigate the general phenomena of spatiotemporal chaos in 
weakly nonlinear oscillatory reaction-diffusion systems. 

To control the CO oxidation reaction, the partial pres- 
sure of CO in the gas above the catalytic surface can be 
modulated in time by varying the dosing rate of CO into 
the reaction chamber. The modulation of the CO pres- 
sure changes reaction conditions in the same way in all 
surface elements and, therefore, its effect is global. By 
making the modulation dependent on the properties of 
monitored patterns, global feedback can be introduced 
in the reaction system. Suppression of chemical turbu- 
lence and development of new, feedback-induced wave 
patterns, including regimes of intermittent spatiotempo- 
ral chaos, have been experimentally observed El 00 an d 
theoretically investigated IEEU fc> r the CO oxidation re- 
action. 

Alternatively, external periodic modulation of the par- 
tial CO pressure can be applied, leading to a periodic 
forcing of this chemical reaction. In experiments [Tlj . 
entrainment of CO oscillations by such periodic forc- 
ing has been detected, but spatially-resolved observations 
of concentration wave patterns on the platinum surface 
were not then possible. Subsequently, spatially-resolved 
experiments on periodic forcing of chemical turbulence 
in this reaction were performed |8|. In addition to full 
entrainment with uniform oscillations, regimes with os- 
cillating cells or labyrinthine patterns, and cascades of 
amplitude defects characteristic of spatiotemporal inter- 
mittency, were observed. Effects of periodic forcing have 
also been experimentally investigated in the oscillatory 
Belousov-Zhabotinsky reaction 0, [ill Q- where uni- 
form oscillations are, however, always stable and sponta- 
neous development of turbulence is not known. 

Previous theoretical studies [HI E3 of fronts sepa- 
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rating different phase-locked states in the periodically 
forced CGLE have shown that, under the conditions of 
the Benjamin- Feir instability, such fronts may undergo 
"explosions" as the forcing amplitude is decreased. When 
such an explosion occurs, the interfacial region separating 
the two locked states grows and eventually fills the en- 
tire medium with a turbulent phase. Similar phenomena 
were also found in coupled map lattice models [13, llSi • 

In this paper, the effects of front explosions under 2:1 
and 3:1 resonant forcing are studied in the realistic model 
of CO oxidation on Pt(llO) for parameter values typical 
for the experiments on this chemical reaction where non- 
harmonic oscillations with substantial amplitudes are ob- 
served. As we show, front explosions occur in very dif- 
ferent ways in these two resonant regimes. Our work 
provides a framework for future experimental studies of 
front explosions in the CO oxidation reaction. 

II. PERIODICALLY FORCED SURFACE 
REACTION 

Model 

We consider the realistic model |3j of catalytic CO ox- 
idation on a Pt(llO) surface. The model takes into ac- 
count adsorption of CO and O2 molecules, reaction, des- 
orption of CO molecules, the structural phase transition 
of the Pt(llO) surface, and the surface diffusion of ad- 
sorbed CO molecules. Letting u(r) and v(r) represent the 
surface coverage of CO and O at position r, respectively, 
and w(r) the local fraction of the surface area found in 
the non-reconstructed lxl form, the reaction-diffusion 
equations giving the time evolution of these fields are 

ut = ki scoPcoO- - u 3 ) - k 2 u - k 3 uv + DV 2 u 
v t = k A po 2 [so, 1x1 w + s , 1x2(1 - w)](l - u - v) 2 
—ksuv) 

All three fields can vary in the interval from [0, 1] and 
we define c(r, t) = (u(r,t),v(r,t),w(r,t)) for future ref- 
erence. For an explanation and specification of the values 
of the parameters see Table UJ Although certain features, 
such as surface roughening, faceting, formation of sub- 
surface oxygen, and the effects of internal gas-phase cou- 
pling are not taken into account, this model has proven 
to be remarkably successful in describing most aspects of 
the experimental observations on this system. In partic- 
ular, it is able to capture oscillations in the CO oxidation 
process seen under some experimental conditions. 

We are interested in situations where the catalytic sur- 
face reaction is subjected to external periodic forcing. 
Experimentally, it is convenient to periodically modulate 
the CO partial pressure pco- Consequently, we assume 
in this paper that pco varies according to the equation 

Pco(t)=Po(l- acosujft), (2) 



Table I: Parameters of the model. 



ki 


3.14 x lO^^mbar" 1 


Impingement rate of CO 


k 2 


10.21s" 1 


CO desorption rate 


k 3 


283.8 s" 1 


Reaction rate 


&4 


5.860 x 10 5 s^mbar" 1 


Impingement rate of O2 


k 5 


1.610 s" 1 


Phase transition rate 


SCO 


1.0 


CO sticking coefficient 


SO, lxl 


0.6 


Oxygen sticking coefficient 






on the lxl phase 


SO, 1x2 


0.4 


Oxygen sticking coefficient 






on the 1x2 phase 


i^o, Su 


0.35, 0.05 


Parameters for the 






structural phase transition 


D 


40/im 2 s" 1 


CO diffusion coefficient 


Po 2 


12.0 x 10" 5 mbar 


O2 partial pressure 


Po 


4.6219548 x 10" 5 mbar 


Base CO partial pressure 



where a is the amplitude of the forcing and uof the fre- 
quency of the forcing. The base value of partial CO pres- 
sure is po. In the present study the parameter values of 
the partial pressures are chosen such that in the absence 
of forcing uniform oscillations are unstable with respect 
to small perturbations and chemical turbulence sponta- 
neously develops in the unperturbed system. The reac- 
tion itself is oscillatory with period To = 2.550049 s for 
a = 0. 

Numerical simulations of the model were performed 
using a first-order finite difference scheme for the spatial 
discretization with a grid resolution of Ax = 4/im. For 
the temporal discretization an explicit Euler scheme with 
a fixed time step At = 0.0001s was used. The system size 
was L 2 = 0.4 x 0.4mm 2 and no-flux boundary conditions 
were imposed. 



Fronts separating resonantly locked states 

We suppose that when periodic forcing is applied to the 
system it is locked in an n : m resonance. In this case 
there are n distinct phases of the oscillation. Chemical 
fronts separating pairs of these n resonantly locked states 
can be formed by a suitable choice of initial conditions. In 
particular, in a two-dimensional domain, we choose initial 
conditions such that the left and the right half planes are 
homogeneous but in two different phases, respectively, 
and separated by a small interfacial zone where the field 
values are chosen at random. 

Under conditions where the unforced system is 
Benjamin- Feir unstable and the forcing amplitude is suf- 
ficiently large, such initial conditions generate a front sep- 
arating the homogeneous phases which has the form of 
an interfacial zone delineated by left and right profiles. 
The situation is depicted in Fig. We let hz J (y,t) = 
min{x : \c(x,y,i) — Cq (t)| > e} and hn(y,t) = max{x : 
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Figure 1: Sketch of a phase front separating homogeneous 
domains of two resonantly locked phases. The interfacial zone 
is delimited by left and right profiles, Jil and Iir, respectively. 



\c(x,y,t) - cft(t)\ > e} (see Fig. [TJ. Here cft(t) and 
Cq (t) are the concentration fields in the two mode locked 
states that the front separates and e is a small numerical 
value which is chosen to be 0.01 in the simulations pre- 
sented below. We can thus define the interfacial zone as 
I(t) = ?/) : h L (y,t) <x< h R (y,t)}. 

The dynamics of the front can be analyzed in terms 
of several quantities that characterize its structure. The 
instantaneous intrinsic thickness of the interfacial zone 
at position y along the front is given by 



A(y,t)=h R (y,t)-h L (y,t). 



(3) 



The spatial average (A)(t) of the intrinsic thickness may 
be computed from 



(A)(4) = 



dyA(y,t), 



(4) 



where, henceforth, the angle brackets will refer to a spa- 
tial average along y. We also define the instantaneous 
position of the front at a point y along the front by 



x f (y,t) = (h R (y,t) + h L (y,t))/2, 



(5) 



and the instantaneous mean position as (xf)(t). After a 
transient time, the phase front dynamics is observed to 
enter a statistically stationary regime where the tempo- 
ral average of the time derivatives of (xf)(t) and (A)(t) 
are independent of time. These quantities are called the 
front velocity Vf and the interface growth rate ^a, re- 
spectively. If va = then the interface has an average 
finite thickness which is denoted by Ao- 

Studies of the 3:1 resonantly forced complex Ginzburg- 
Landau equation have shown that if the amplitude of 
the forcing is reduced, a critical value is reached where 
the average thickness of the interfacial zone Ao is no 
longer finite. For such values of the forcing, when the 
system starts from the two-phase initial conditions de- 
scribed above, the interfacial zone grows to fill the entire 



domain. We term this phenomenon a "front explosion" . 
We now investigate this phenomenon for the resonantly 
forced catalytic surface oxidation reaction. 



III. 2:1 FORCING 

In this section we study the front dynamics when the 
system is close to the 2:1 resonance, i.e., Tf = To/2 and, 
hence, ujf = 2coo. In the 2:1 resonance regime we have 
two distinct equivalent phases that differ by a phase shift 
of 7r. An example of a front separating the two phase 
states is shown in the left panel of Fig. [3 The front is 




Figure 2: Snapshots of u(r) in the asymptotic state for Tf = 
1.275s without front explosion (left, a = 0.0138) and with 
front explosion (right, a = 0.00864). Black corresponds to 
small values of u and white to large values of u. 

stationary with finite average width Ao. As the ampli- 
tude of the forcing is decreased we observe that Ao grows 
and, for amplitudes lower than a critical value, the inter- 
facial zone grows without bound and the front explodes. 
An example of an exploding front is shown in Fig. 03 In 




Figure 3: Snapshot series of u(r) for a — 0.0130 showing a 
typical front explosion. Left, t = 5 X T/; right, t = 20 x Tf 
with Tf = 1.275s. Color coding as in Fig. El 

this case the turbulent zone grows to fill the finite do- 
main leading to turbulent dynamics in the entire system 
as shown in the right panel of Fig. 

The turbulent state that emerges when the interfacial 
zone explodes has the characteristics of defect-mediated 
turbulence, where the dynamics of a pattern is dominated 
by the rapid motion, nucleation and annihilation of point 
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defects (vortices or dislocations) [Oj . A defect is charac- 
terized by its integer topological charge (or winding num- 
ber) m top which is defined by ^ § V0(r, i) • dl = ±m top 
[20j | . where 0(r, t) is the local phase and the integral is 
taken along a closed curve surrounding the defect. A 
topological defect, thus, corresponds to a point in the 
medium where the local amplitude is zero and the phase 
is not defined. In the CGLE only topological defects 



with rn to 



±1 are observed and that is also the case 



here. In Fig. 0] (left panel) we see that the number of 
defects fluctuates about a statistically stationary aver- 
age value. The probability distribution of defects with 
positive topological charge shown in the right panel of 
this figure has an approximate Gaussian form. This is 
in accord with a simple model of defect dynamics based 
on rate equations and taking into account no-flux bound- 
ary conditions 0, [2^ . The model predicts a modified 
Poissonian distribution which, in the limit of large mean 
defect numbers, converges to a Gaussian. 
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Figure 4: (Color online) Left: Number of defects n to tai as 
a function of time (same forcing as in the right panel of 
Fig. for 0(r, t) = arctan[(w(r, t) - w°)/(u(r, t) - u )] with 
(w°,u°) = (0.354,0.487) as the center of rotation. Right: 
Probability distribution of the number of defects with posi- 
tive topological charge n+ (red circles). It is clearly different 
from a squared Poissonian (black triangles) |2l| but shows 
very good agreement with a Gaussian (blue diamonds). 

A quantitative characterization of the transition from 
a stationary front with finite width (vf = va = 0) to an 
exploding front (vf = and va > 0) is given in Fig. [51 
which plots both Ao and va versus the forcing ampli- 
tude a. The transition occurs at a critical value a\ which 
coincides, within the numerical resolution, with a bifur- 
cation leading to the coexistence and bistability of ho- 
mogeneous oscillatory states and turbulent states. The 
homogeneous oscillatory states become unstable only for 
values of a much smaller than a^, namely at a ~ 0.0006. 

There are number of interesting features of the dy- 
namics as the system approaches the front explosion 
point. Approaching a\ from above, the front explo- 
sion is preceded by a period-doubling bifurcation at a^. 
Fig. shows that the bifurcation occurs in the homo- 
geneously oscillating domains as well as in the interfa- 
cial zone. This is in accord with the behavior of a sin- 
gle oscillator close to the 2:1 resonance. The phase dia- 
gram showing its states as a function of the forcing fre- 
quency and amplitude near the 2:1 resonance is given 
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Figure 5: (Color online) Plot of the interface width Ao (red 
circles) and the interface growth rate va (blue stars) as a 
function of the forcing amplitude a for Tf = 1.275s. The 
statistical errors are less than the size of the symbols. Three 
different regimes can be identified. See text for details. 



in Figure As the forcing amplitude is increased at 
oof/coo = 2 one crosses a parameter region where period 
doubling occurs. For the spatially distributed system, 
within this period-doubled or 4:1 regime, the interface 
dynamics (c)j = Aj 1 J /( ^ dr c(r, t) undergoes further 
period-doubling bifurcations into chaos before the front 
explodes (see Fig.EJ. Here A/ is the area of the interfacial 
zone. Note that in the 4:1 regime two different types of 
fronts can exist, namely tt and tt/2 fronts. Here, we have 
focused on the tt/2 fronts since they directly correspond 
to the (tt) fronts present in the 2:1 regime. Simulations 
show that 7r fronts do not explode at a = a\. They are 
stationary even below a\ and do not become unstable 
within the 4:1 regime. 

Phase fronts can also be observed for random initial 
conditions. As long as a > aj, the turbulent phase is 
unstable and homogeneous domains separated by phase 
fronts are nucleated. If these fronts are sufficiently curved 
they are not stationary but propagate. In a finite sys- 
tem as considered here, this can eventually lead to the 
generation of a spatially homogeneous state without any 
fronts. Interestingly, tt/2 fronts survive much longer than 
tt fronts. 

Most of the results presented in this section do not 
qualitatively depend on the exact forcing frequency as 
long as ujf is close to 2coo. Restricting ourselves to the 
range of forcing frequencies present in Fig. [7[ only the 
features of the dynamics as the system approaches the 
front explosion point from above depend on cof. If ujf 
is too small, the period doubling of the dynamics within 
homogeneous domains and interfacial zones is absent. If 
ou f is too large, a cascade of period doubling bifurcations 
can be observed in the homogeneous domains as well. In 
general, a\ increases with increasing detuning from the 
resonance. 
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Figure 6: (Color online) Dynamics of the fields c(r) = 
(u(r),v(r),w(r)) averaged over the interfacial zone /(£), de- 
noted as (c)j, for Tf = 1.275 and different values of the forc- 
ing amplitude a > corresponding to different states of the 
system (a = 0.0173 (2:1), a = 0.0155 (4:1), a = 0.0147 (8:1), 
a = 0.0138 (chaos)). The interface dynamics undergoes a 
period-doubling cascade into chaos while the dynamics of the 
homogeneous oscillations outside the interfacial zone (solid 
black curve) has only a single period-doubling bifurcation at 
a% ~ 0.0162. The limit cycle of the homogeneous oscillations 
is basically the same for the subngures denoted 4:1, 8:1 and 
chaos. 



IV. 3:1 FORCING 

Near the 3:1 resonance where Tf = To/3 (ujf = 3coo) 
the front dynamics exhibits a number of features which 
differ from those described above for the 2:1 resonance. 
In the 3:1 resonance regime the broken translational sym- 
metry gives rise to three distinct phases that differ by 
a phase shifts of 27r/3; consequently, fronts separating 
any two of these three phases will propagate with fi- 
nite velocity Vf. For sufficiently large forcing ampli- 
tudes, these propagating fronts are nearly planar with 
only small transverse structure. 

As the forcing amplitude decreases, the front width 
Ao increases in magnitude and, for amplitudes lower 
than a critical value a%, the interfacial zone grows with- 
out bound and the front explodes. When the interfacial 
zone exhibits unbounded growth, the mean velocity of the 
front is zero, Vf =0. A simulation of an exploding front 
is shown in top two panels of Fig.[SJ In contrast to the ex- 
ploding front in the 2:1 resonantly forced regime, the ex- 
panding interfacial zone does not display spatio-temporal 
turbulent dynamics, instead, a labyrinthine patterns de- 
velops. This can be seen in bottom two panels of Fig. [S| 
which shows the formation of a labyrinthine pattern after 
the interfacial zone fills the simulation domain. Examina- 
tion of the local dynamics indicates that the labyrinthine 
pattern is locked 2:1 to the applied forcing. 

Insight into origin of this behavior can be obtained 




Figure 7: (Color online) 2:1 Arnold tongue for a single oscilla- 
tor. Note the appearance of a period-doubling cascade within 
the tongue — well below the 1:1 tongue. For cof/ooo = 2, 
only the first period-doubling bifurcation is present. Its lo- 
cation compares very well with the observed location for the 
spatially extended system. 



by considering the dynamics of a resonantly forced sin- 
gle oscillator. The phase diagram showing its states as 
a function of the forcing frequency and amplitude near 
the 3:1 resonance is given in Fig.JUJ This phase diagram 
has a number of features that differ from those of the 
2:1 Arnold tongue shown earlier in Fig. [7| The period 
doubling region within the Arnold tongue does not in- 
tersect the uo f = 3ujq line. In addition, there is a region, 
indicated by the cross pattern, where there is bistability 
between 3:1 and 2:1 locking. In the spatially extended 
system, the values of the forcing amplitude used to obtain 
the exploding front are well within the area of bistability 
and coexistence of 3:1 and 2:1 states for a single oscil- 
lator, suggesting the possibility of appearance of a 2:1 
resonantly locked state as observed in the simulation. Of 
course, the single oscillator phase diagram does not pro- 
vide information on the spatial structure of the pattern 
and the Arnold tongue for the spatially distributed sys- 
tem may display a richer structure. Our simulation result 
confirms the existence of a 2:1 state with a labyrinthine 
pattern within the 3:1 resonance tongue. 

The front explosion can again be characterized quan- 
titatively by plotting the mean interfacial width, Ao, the 
growth rate of the mean interfacial width, v&, and the 
mean front velocity, Vf as a function of the forcing ampli- 
tude. The data in Fig. ITU1 shows that there is a transition 
from a propagating front with Vf ^ and va = to an 
exploding front with Vf = and va > at a critical 
value a^. 

This figure also shows that there are a number of in- 
teresting features of the dynamics in the vicinity of the 
front explosion. Approaching a\ from above, the front 
explosion is preceded by a transition of the front itself. 
At af, the front velocity suddenly drops while the width 
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Figure 10: (Color online) Plot of the mean interface width Ao 
(open circles), the mean front velocity Vf (open diamonds) 
and the interface growth rate va (stars) as a function of the 
amplitude of the forcing a for Tf = 0.85s. The statistical 
errors are less than the size of the symbols. Three different 
regimes can be identified. See text for details. 



Figure 8: Series of snapshots of the u(r, t) field for a = 0.0354 
and Tf = 0.85s showing the exploding front and the eventual 
formation of a labyrinthine pattern. From left to right and 
top to bottom: t = 40 X T f , t = 760 x 7>, t = 4360 x T f , 
t = 10600 x Tf. Color coding as in Fig. Note that due 
to the initial conditions and the no-flux boundary conditions 
certain symmetries persist in the pattern. 
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Figure 9: (Color online) 3:1 Arnold tongue for a resonantly 
forced single oscillator. The period-doubling bifurcation 
within the tongue is shown as a region shaded with paral- 
lel lines. Note that this region is does not intersect the 3:1 
resonance line (vertical line in the figure) in contrast to the 
case of 2:1 forcing. The transition from 3:1 to 2:1 locking 
(solid line with *) is subcritical and the region of bistability 
is marked by the orange cross pattern. 



of the inter facial zone increases. As a is reduced for 
CL3 < a < G&3, the width increases and the front veloc- 
ity decreases up to the transition point a\ where Vf = 0. 
The sudden change in Vf and Ao at « 0.038546 is 
signalled by behavior of the transient time needed to es- 
tablish the interface. Figure ITT1 plots the temporal evolu- 
tion of the average width and location of the interfacial 
zone for Tf = 0.85 and different values of the forcing 
amplitude. For both values of the forcing amplitude two 
regimes may be distinguished in the figure. For short 
times a relatively thick and slowly propagating front is 
observed which persists for a transient time t c indicated 
by the dashed lines in the figure. For times t > t c , the 
front transforms into a faster moving thin front. The 
evolution of the interfacial zone for t < t c is very similar 
to that observed for a < a$. As a approaches from 
above, the transient time t c increases as shown in the 
figure. For a < af, no evidence for a finite t c has been 
found. Thus, these data suggest that transient time t c 
diverges at leading to a distinctly different front types 
below and above a%. 

Another interesting feature of the dynamics in the 
vicinity of the front explosion is that, approaching a\ 
from above, the front explosion is preceded by a sequence 
of period-doubling and inverse period-doubling bifurca- 
tions in the dynamics of the interfacial zone characterized 
by the behavior of (c)/(t). Depending on the value of a, 
(c)j(t) is locked either 3:1 or 6:1 to the forcing; however, 
there is no period-doubling bifurcation of the homoge- 
neous oscillations outside the interfacial zone in contrast 
to the case of 2:1 forcing. This observation is expected, 
based on the phase diagram for a single oscillator given in 
Fig. [HI which does not show a period-doubling bifurcation 
close to the 3:1 resonance. 

Finally, we observe that if random initial conditions 
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Figure 11: (Color online) Temporal evolution of the average 
width (A) and position (xf) of the interfacial zone for Tf = 
0.85 and different values of the forcing amplitude a > af. The 
thin light (green) curves correspond to a = 0.03887 and the 
thick black ones to a = 0.03865. Two regimes are visible. At 
short times, a thick and slowly propagating front is present 
which transforms into a thin and fast propagating front at 
t c . The transient time t c is shown as the dashed lines and 
increases with decreasing a. 

are chosen, labyrinthine patterns form for a < while 
homogeneous domains locked 3:1 to the forcing and sep- 
arated by phase fronts occur for a > a%. In the latter 
case, due to the propagation of the fronts, these inhomo- 
geneous states generally evolve to homogeneous asymp- 
totic states in a finite system. 

V. CONCLUSIONS AND DISCUSSION 

The results presented in this paper show that while 
front "explosions" occur in a realistic model for catalytic 
CO surface oxidation reactions, a number of new phe- 
nomena are observed. In the 3:1 regime of the resonantly 
forced CGLE, when the underlying unforced system is 
Benjamin- Feir unstable [TH L fT ^I . and also for the period-3 
piecewise linear coupled maps |l7l , propagating fronts 
have a turbulent interfacial zone separating homogeneous 
resonantly- locked states. The interfacial zone grows in 
this case as a power law when the magnitude of the forc- 
ing amplitude is decreased. In the CO oxidation reaction, 
on the other hand, the interfacial zone in the propagat- 



ing front is much more structured and front explosion 
leads to the formation of a labyrinthine pattern with 2:1 
locking. Investigations of the front explosion dynamics 
in the 2:1 regime of the CGLE with the Benjamin- Feir 
instability reveal propagating Bloch fronts with a tur- 
bulent interfacial zone which grows without bound be- 
yond a certain critical value of the forcing amplitude [2^ . 
In both the 2:1 CGLE and in period-2 coupled maps, 
poorly-characterized front explosions of stationary fronts 
have also been observed. In the 2:1 resonantly forced 
CO surface oxidation reaction, the stationary front sepa- 
rating the two homogeneous resonantly locked states has 
a turbulent interfacial zone. As the forcing amplitude 
is decreased, beyond a critical value the interfacial zone 
grows without bound leading to a defect-mediated tur- 
bulent state. 

Although spatially-resolved experiments with periodic 
forcing of the CO oxidation reaction have already been 
performed the behavior of fronts was not studied. 
To create a front, special initial conditions are needed, 
which were not implemented in these experiments where 
forcing was always applied starting from the turbulent 
initial state. An additional difficulty was that, in the 
experiments, the frequency of uniform oscillations in ab- 
sence of forcing could not be exactly measured, because 
such oscillations were unstable and led to turbulence. In 
these experimental investigations, the attention was fo- 
cused on qualitative properties of different patterns in- 
duced by forcing. Remarkably, irregular stripe patterns 
with the 2:1 locking were found then for high forcing 
frequencies corresponding to 3:1 forcing |8|. Such ir- 
regular stripe patterns were interpreted as labyrinthine 
patterns affected by strong surface anisotropy. In the 
present theoretical study, we have indeed found that the 
2:1 locked labyrinthine patterns should exist under 3:1 
resonant forcing in the CO oxidation system. 

Our analysis of front explosions provides a theoreti- 
cal framework for future experimental studies of the CO 
oxidation reaction under resonant periodic forcing. Gen- 
erally, it provides an example of the front explosion phe- 
nomena for a system with non-harmonic oscillations and 
more complex dynamics, different from the ideal situa- 
tion described by the complex Ginzburg-Landau equa- 
tion. 
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